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Abstract 

We present a measurement of the Kl lifetime using the KLOE detector. From a 
sample of ~ 4 x 10^ KsK^ pairs following the reaction e~^e~ ^ — > KsKl we 
select ~ 15 X 10^ tt^tt^tt^ decays tagged by Ks — > tt^tt" events. From a fit 

of the proper time distribution we find tl = (50.92 it 0.17stat i 0.25syst) ns. This is 
the most precise measurement of the Kl lifetime performed to date. 

PACS: 

keywords: 



1 Introduction 



The Ki lifetime is necessary to determine its semileptonic partial widths from 
the branching ratios (BR). The partial widths can be used to extract the 
CKM matrix element \Vus\- Present knowledge of t{Kl) comes from a single 
measurement performed more than 30 years ago [1] and its error dominates 
the uncertainty in the partial Kl decay rates. At DA$NE , the Frascati 0- 
factory , nearly monochromatic i^'j^-mesons are produced with p~110 MeV/c 
corresponding to a mean path of 340 cm. The KLOE detector is large enough, 
r=200 cm, so that ~50% of the Kl decay inside it. The statistical error on the 
lifetime depends strongly on the time interval covered in the measurement [2]: 
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where T = At/r is the time interval observed, in KL-\iietime units. With 
T~0.4 and A^~9 x 10^ we can reach an accuracy of ~0.3%. 

We have measured the hfetime using the decay ti^ti^ti^ tagged by 

Ks-^TT'^iT^ events. This choice maximizes the number of usable events and 
minimize the disturbance of the decay on the detection of the tagging Ks 
decay and therefore the systematic uncertainty. 



2 Experimental setup 

In DA$NE electrons and positrons collide with an angle of 25 mrad and a 
center of mass (CM) energy W = M{(f)). 0-mesons are produced with a cross 
section of ~3 fib and a transverse momentum of ~ 12.5 MeV/c toward the 
center of the collider rings. The energy W, the position of the beam crossing 
point {x, y, z) and the momentum are determined from Bhabha scattering 
events. In a typical run of integrated luminosity / £dt~100 nb~\ lasting about 
30 minutes, the corresponding errors are: 5W=AQ keV, 5p<^=30 keV/c, 5x = 
30 fim, and 6y = 30 fim. 

The detector consists of a large cylindrical drift chamber, DC [3], whose axis, 
defined as the z-axis, coincides with the bisectrix of the two beams. The DC 
is surrounded by a lead-scintillating fiber sampling calorimeter, EMC [4] . The 
DC and EMC are immersed in a solenoidal magnetic field of 0.52 T with the 
axis parallel to the beams' bisectrix. The DC tracking volume extends from 
28.5 to 190.5 cm in radius and is 340 cm long. The transverse momentum reso- 
lution is (5p_l/p_l~0.4%. Vertices are reconstructed with a resolution of ~3 mm. 
The calorimeter is divided into a barrel and two endcaps and covers 98% of the 
solid angle. Photon energies and arrival times are measured with resolutions 
as/E = 0m7/^E (GeV) and at = 54 ps/^E (GeV) © 50 ps respectively. 

Photon entry points are determined with an accuracy cr^ ~ 1 cm/ E (GeV) 
along the fibers and (T_l~1 cm in the transverse direction. A photon is defined 
as an EMC cluster of energy deposits not associated to a track. We require that 
the distance between the cluster centroid and the entry point of the nearest 
extrapolated track be greater than 3cr, a = az (B (J±. 

The trigger [5] uses information from the calorimeter and chamber. The EMC 
trigger requires two local energy deposits above threshold {E > 50 MeV in the 
barrel, E > 150 MeV in the endcaps). Rejection of cosmic-ray events is also 
performed at trigger level, checking for the presence of two energy deposits 
above 30 MeV in the outermost calorimeter planes. The DC trigger is based 
on the multiplicity and topology of the hits in the drift cells. The trigger has 
a large time spread with respect to the beam crossing time. It is therefore re- 
synchronized with the machine radio frequency divided by four, Tgync = 10.85 
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ns, with an accuracy of 50 ps. During the 2001-2002 data taking the bunch 
crossing period at DA$NE was T=5.43 ns. The correct colhsion time, To, of 
the event is determined off-hne during event reconstruction [6]. 



3 Data analysis 

</)-mesons decay into Ks-Kl pairs ~34% of the time. Production of a Kl 
is tagged by the observation of i^5— ^■7r"'"7r~ decay. The Kl—>^ decay 
vertex is reconstructed along the direction opposite to that of the Ks in the (p 
rest frame. The chamber alone measures the Ks^ti^ti^ decay and therefore 
the direction of the K^. The decay vertex and the photon energies are 
obtained from EMC information. The data sample, collected during 2001 and 
2002, corresponds to an integrated luminosity of ~400 pb~^. Some 1.2 x 10^ 
0-mesons were produced. Additional details can be found in reference [2]. 
Ks — TT^TT" decay events must satisfy the following requirements: 

(1) There must be two tracks with opposite charge, forming a vertex V in 
a cylinder with ry < 10 cm, \zy\ < 20 cm. No other tracks should be 
connected to the vertex. 

(2) The Ks momentum in the rest system, must satisfy 100 < pKg < 120 
MeV/c. The 7r"'"7r~ invariant mass M^tttt) must satisfy 492< M^tttt) <503 
MeV/c2. 

The efficiency for finding Ks tt'^tt^ events is e ~ 68%. The position of the 
(j) production point, x^, is determined as the point of closest approach of the 
Ks momentum, propagated backwards from the Ks vertex, to the beam line. 
The Ks — >■ vr+TT" decay provides an almost unbiased tag for the K^ when it 
decays into neutral particles and a good measurement of the K^ momentum, 
PKj^ = Pit) — PKs- The accuracy in the determination of the Kl direction is 
obtained from Kl^-k^tx^-k^ events, measuring the angle between px^^ and the 
line joining the (f) production point and the tt^tt'tt^ decay vertex. We find 
= 1.5°, ae = 1.8°. 

The position of the Kl vertex for Kl^t^^tx^tx^ decays is obtained from the 
photon arrival times at the EMC. Each photon defines a triangle CDE, see 
figure 1 (left), where Ik is the Kl path length, is the distance from the Kl 
decay point D to the entry point E and d is the distance from the cluster to 
the collision point C. From the known positions of C and E, the ECD = 9 
angle and the time spent by the kaon and the photon to cover the path CDE 
we find the length of CD. There are two solutions. One has D along the Ks 
path and is rejected. The position D of the Kl decay vertex is obtained from 
the energy weighted average of the two closest lK,i, Ik = Hi}K,i x Ei)/J2Ei 
where i is the photon index. Finally we require at least one third photon with 
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Ik,3 -Ik\ < 5 X a{lK)- 




Fig. 1. Left. The CDE triangle. Right. Distribution of the difference iRi"^^) 
— /i^(7r^7r~) for — > 7r+7r~7r'^ events as a function of /^^-(Tr+vr"). See text. 

The accuracy of the Ik determination is checked by comparing the path 
measured by timing with the calorimeter and, with a much better accuracy, 
by tracking with the DC, for i<'i^7r"*"7r~7r° decays. The path length from the 
calorimeter timing has on average a constant offset of 2 mm with respect 
to the value obtained with the DC, figure 1, right. The determination of Ik 
depends crucially on the correct identification of the collision time Tq. Using 
again i^^— i>7r"'"7r~7r° events we have verified that Tq is incorrect less than 0.1% 
of the time. 

The resolution (j{Ik) is determined from 7r"'"7r~7r° events by comparing 

^xiT^^) and /x(vr"''7r~), where the former is the weighted average obtained from 
the two photons from vr" and (vr"''-/?") is the distance between the vertex 
of the two charged decay pions and the production point. An example of 
the /A'(7r°) — Ik^t^'^it') distribution is shown in figure 2, left. It has been 
fitted both with a single and a double gaussian [2]. In the case of the double- 
gaussian fit the relative weights of the two components are free parameters 
of the fit. The single-gaussian fit gives an average resolution of ~ 2.5 cm. In 
the double-gaussian fit, the bulk of the distribution (~ 82% of the total) has 
0"! ~ 2.1 cm, while the broader part (~ 18%) is well described by a gaussian 
with (72 ~ 5.4 cm. The behaviour of the resolution has been studied as a 
function of /i^(7r"^7r~): we find a quadratic dependence on Ik in both cases 
(single and double-gaussian fit). For the single-gaussian fit we have (t{Ik) = 
1.65 + 0.59 X 10-2 xIk + 0.45 x 10"^ x Ij^ cm (figure 2, right). In the double- 
gaussian fit the relative weights of the two components change as a function 
of Zx(7r+7r-). Since for each point the population weighted average of cxi and 
(J2 agrees at the 10% level with the a of a single-gaussian fit, we use the latter 
as estimate of the vertex resolution. The effects of the tails on the fit value 
are discussed in Section 6. 



5 




-30 



-20 



-10 10 20 30 




20 40 60 80 100 120 140 160 180 

l„(7c'^7t"\cm) 



Fig. 2. Left. Distribution of the difference Ik{t^^) — Ik{'^~^'^~) for tt^tt^tt^ 
events. Right. The a obtained with a single-gaussian fit as a function of Ixi'^'^'^")- 
See text. 

The tagging efficiency has been evaluated by MC as a function of Ik for 
the dominant Kl decay channels. The results are shown in figure 3. The 
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Fig. 3. Tagging efficiency as a function of Ik for the main Kl decay modes. 

difference in tagging efficiency among the Kl decay modes is mainly due to 
the dependence of the trigger efficiency. Only the calorimeter trigger [5] is used 
for the present analysis. The trigger efficiency is, on average, ~100% for Kl —>■ 
^o^o^o between 85 — 95% for charged Kl decays. The trigger efficiency 
also depends on the position of the Kl decay vertex. Another contribution is 
the dependence of the reconstruction efficiency for the pion tracks from Ks — > 
7r"'"7r~ on the presence of other tracks in the drift chamber. This contribution 
depends on the position of the Kl decay point and affects mainly events with 
Kl — * charged particles near the production point. 



The tagging efficiency for the Kl 
pendence on Ik, with a slope of b 
a = (68.04 ±0.01)%. 



TT^TT^Tr^ channel has a small linear de- 
= (1.2 ± 0.2) X 10~^/cm, and a constant 
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4 Kl 7r''7r"7r'' acceptance 



The KL-^TT^TT^ir^ decay has a relatively large BR, ~21%, and has very low 
background. KL-^TT^TT^'ir^ events are accepted if at least three calorimeter clus- 
ters are found satisfying: 

(1) Energy larger than 20 MeV. 

(2) Distance from any other cluster larger than 50 cm. 

(3) No association to a chamber track. 

(4) j — Ik\ < 5 X a^lx), where Ik is the energy weighted average of the 
two values of lK,i nearest together. 

For the Kl lifetime measurement, we retain events with 40 < Ik < 165 cm 
and a polar angle 6 in the interval {40°, 140°}. These two conditions define 
the fiducial volume (FV). The main sources of event losses are: 1) geometrical 
acceptance; 2) cluster energy threshold; 3) merging of clusters; 4) accidental 
association to a charged track; 5) Dalitz decay of one or more 7r*^'s. The effect 
of these inefficiencies is to modify the relative population for events with 3, 4, 
5, 6, 7 and > 8, clusters with a loss of global efficiency of ~ 0.8%. 

Monte Carlo (MC) simulations, based on the KLOE standard MC [6] show 
that event acceptance with the above selection has a linear dependence on Ik, 
e{lK) = (0.9921 ±0.002)- (1.9 ± 0.2) X 10"^ x/^ with Ik in cm, (figure 4, left) 
mainly due to the vertex reconstruction efficiency. This has also been checked 
using Kl —>■ n~^7r~n^ events both from data and MC. We find the same lin- 
ear dependence, with compatible slopes within their statistical uncertainties, 
figure 4, right. 





102 


g 








>> 
u 


101 


a 








fici 


100 












99 




98 




97 




96 





'ficiency (% 




M 96 


\ 


94 


1 


92 
90 


i 

175 200 


88 



Ik i^^) 



25 50 75 100 125 150 175 200 

1^ (cm) 



Fig. 4. Left: Vertex reconstruction efficiency as a function of the decay path length 
for Kl — > 7r'^7r''7r'^ Monte Carlo events. Right: the same for Kl — > 7r"''7r~7r'' data 
(triangles) and Monte Carlo (squares) events. 



A comparison between data and MC of the photon multiplicity and total en- 
ergy distributions for Kl — 7r°7r°7r° decays shows that only events with three 
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and four clusters contain some background. Background, mostly to the three 
cluster events, is due to — > ti^tt~tt^ decays where one or two charged pions 
produce a cluster not associated to a track and neither track is associated to 
the Kl vertex. Other sources of background are ti^tt^ decays (possibly 

in coincidence with machine background showering close to the collision point 
generating soft neutral particles) and Ks tt^tt^ following Kl — > Ks regen- 
eration in the DC material. The Kl tt^tt^tt^ background and the other 
backgrounds in the three cluster sample are strongly reduced by requiring at 
least one cluster in the barrel with E >50 MeV and no tracks approaching the 
Kl line of flight by less than 20 cm. The efficiency vs Ik for the three cluster 
sample has been found by MC. It is almost flat with an average ~55% inside 
the FV. 

In Table 1 we show the fractions of the background components before and 
after the background cuts. The background contamination is reduced from ~ 
4.9% to ~ 1.3% with an efficiency on the signal of ~ 99.6%. The distributions 
of the total photon energy for events with 3, 4, 5, 6, 7, > 8 photons are shown in 
figure 5. For three and four photon-cluster samples the different contributions 
from the residual background components are also shown. 

The fractions of events with = 3, 4, 5, 6, 7, > 8 in the FV are given in Table 
2, together with MC results. The few percent differences between data and 
MC are mostly due to a higher proportion of split clusters in the MC than in 
the data. 





N-y >3 before cuts 


Ny >3 after cuts 


Kl channel 


Events 


B/(S+B) 


Events 


B/(S+B) 


signal+backgrounds 


10,536,674 




10,114,899 




all backgrounds 


518,520 


4.92 % 


133,535 


1.32 % 


Kl 'K~^1T~'K^ 


325,076 


3.08 % 


44,917 


0.44 % 


Kl 'nfiv 


28,917 


0.28 % 


3,583 


0.03 % 


Kl — > TTev 


49,140 


0.47 % 


6,062 


0.06 % 


Kl ^ vrO^o 


43,436 


0.41 % 


42,313 


0.42 % 


Kl^Ks^ 7r°7r° 


30,273 


0.29 % 


28,166 


0.28 % 


Kl other 


41,298 


0.39 % 


8440 


0.08 % 



Table 1 

Event types in the fiducial volume from Monte Carlo before and after background 
cuts. The background contamination B /{S + B) is also shown. 
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Total Energy for N=6 (MeV) Total Energy for N=7 (MeV) Total Energy for N>8 (MeV) 



Fig. 5. Kl vr vr vr selection: distribution of the total energy for events with 3, 
4, 5, 6, 7 and > 8 photon clusters. Dots are data, solid histogram is Monte Carlo 
simulation for Kl all channels. Monte Carlo histograms are normalized to the 
same number of entries for data. 



Number of clusters 


Data 


Monte Carlo 


3 


1.163 ± 0.004 % 


0.980± 0.003 % 


4 


7.64 ± 0.01 % 


7.01 ± 0.01 % 


5 


30.22 ± 0.02 % 


28.65± 0.02 % 


6 


57.77 ± 0.03 % 


60.12± 0.03 % 


7 


3.091 ± 0.006 % 


3.074± 0.001 % 


> 8 


0.106 ± 0.001 % 


0.151± 0.001 % 



Table 2 

Fraction of events with 3,4,5,6,7 and > 8 neutral clusters connected to the Kl 
decay vertex in data and Monte Carlo. 

5 Fit of the proper time distribution 



The Kl proper time, t*, is obtained event by event dividing the decay length 
by P'j of the in the laboratory, t* = Ik/{P'Jc). In figure 6 we show the 
t* distribution obtained with ~ 14.7 x 10^ tagged Kl — > events. The 

residual ~ 1.3% background is subtracted using MC results. The variation of 
the vertex reconstruction efficiency as a function of the decay length is taken 
into account by correcting bin by bin the t* distribution with product of the 
MC Kl — > Svr'' efficiency (figure 4, left) and the data/MC efficiency ratios for 
Kl 7r+7r~7r° (figure 4, right). The statistical uncertainty (~ 0.1%) of the 
efficiency estimate is included in the error. 
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Both the background subtraction and the vertex reconstruction efficiency cor- 
rection affect the number of events per bin at the ~ 1% level and the combined 
effect of both corrections leaves the effective statistics essentially unchanged. 
Figure 6 is therefore representative of the sample statistics. 

The t* distribution is fitted with an exponential function over the range 6< 
t* <24.8 ns. This corresponds to a time interval T = At*/r~0.37. With 
~ 8.5 X 10^ events in the fit region we obtain: 

T = (50.87 ±0.17) ns 

with a x^- value of 58 for 62 degrees of freedom (figure 6). 



"0 5 10 15 20 25 30 ^ 35 5 10 15 20 25 30 

t* (ns) t" (ns) 

Fig. 6. Fit of the proper time distribution (left) and residuals of the fit (right). 
Crosses are data and solid histogram is Monte Carlo. The fit is shown as the thick 
sohd line. 

6 Systematic uncertainties 

The number of Stt*^ decays at the end of the selection is given by: 

iV3.o(/i^) = iV3,o(0) J etotil'K) X e-''Kl>^' xgilK - I'k)^'^, + N,,k{lK), (2) 

where Nbck is the residual background at the end of the signal selection, etotil'K) 
is the signal efficiency (etot = etag x Csei) and gilx — I'k) is the vertex resolution 
function. Finally A^; = (l/A^ + l/A/)~^ is the effective mean decay length 
taking into account the interactions inside the chamber (in the gas mixture 
and wires) and A^ is the mean decay length. 

Many effects distort the proper time distribution and have been corrected 
for. The uncertainty in the corrections is included in the systematic error on 
the Kl lifetime. As noted in Section 3, the tagging efficiency (figure 3) is 
well described by a linear function of Ik with a constant term a = (68.04 ± 
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0.03)% and a slope b={1.2 ± 0.2) x 10~^/cm. While the tagging efficiency 
is easily parametrised, it has a significant effect on the overall statistics of 
the sample. Therefore, the value of the lifetime is corrected for the effects 
of the Ix dependence of the tagging efficiency using an analytical correction: 
A* ~ A(l + (h/a)\) [2]. This results in a correction on the lifetime of —0.6% 
with a systematic uncertainty of ±0.1%. 

We also vary the threshold of the cluster energy of the pions from Ks from 
40 to 70 MeV in 10 MeV steps. The slope changes by 0.5% with a systematic 
uncertainty of ±0.25%. 

As discussed, the vertex reconstruction efficiency has been corrected for its 
dependence on l^- We assign a systematic error of ±0.2% due to the statistical 
uncertainty on the slope of the data/Monte Carlo efficiency ratios evaluated 
with Kl — s> 7r'^7r~7r° events. 

We investigate the effects of the cluster energy threshold for photons Ethr, 
by varying Ethr from 10 to 35 MeV in steps of 5 MeV and repeating the full 
analysis. The value of this threshold affects dramatically both the background 
contamination and the relative weights of the samples of different photon- 
cluster multeplicity. For example, if Ethr goes from 20 MeV to 15 MeV, the 
relative weight of the three photon-cluster sample is reduced by almost a factor 
2 while the background increases by 20 %, affecting mainly the four photon- 
cluster sample. Nevertheless, the fit changes by < ±0.2% for 15 < Ethr < 
35 MeV, which we take as a systematic uncertainty. The total systematic 
uncertainty due to the event selection is therefore ±0.3%. 

The effect of the vertex resolution on the fit value has been studied by smearing 
the values sampled from an exponential function with the measured ex's (as 
a function of Ik) both in the case of a single-gaussian fit and in the case of 
a double-gaussian fit [2]. In the second case, the smearing is performed by 
taking into account the relative weights of the two gaussians at a given Ik- 
When the fit is performed to a generated sample of Ik values without or with a 
smearing with the known la resolution parameters, the lifetime value changes 
by well under 0.1 %. Resolution effects are thus negligible in determining the 
final value of lifetime. However, if the resolution parameters a are ~ 10% 
larger than the measured values, the effect on the value of the Kl lifetime is 
~ 0.1%. The (j's are know at 1-2 % level. For the systematic uncertainty in the 
lifetime value for vertex resolution effects we assign a symmetric 0.1 % error, 
which is conservately based on the assumption that the resolution parameters 
can be underestimated by as much as 10 %. 

Kl interactions with the material inside the chamber bias the lifetime mea- 
surement since they reduce the Kl mean path by (1 ± A^/A/)"^ ~ (1 — \l/^i) 
[2]. The interaction rates for regeneration and A or E production are deter- 
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mined from data [7] . The contribution of the regeneration in the DC material 
is found to be ~x3 times lower in data than in the MC prediction. In data 
the contribution of the total nuclear interactions is ~ 0.33% to which we as- 
sign a conservative error of 50%, (0.33±0.16)%. Therefore, the lifetime is 
corrected by +0.33% with a systematic uncertainty of 0.16%. 

The background (table 1) is subtracted from the proper time distribution using 
MC simulation. A residual correction of +0.2% due to background subtraction 
has been added due to the fact that the background from the Ks regeneration 
in the drift chamber material is a factor of three higher in the Monte Carlo 
than in data [7]. Therefore, from Table 1, the 0.28% contribution from the 
reaction must be reduced to 0.1% and the global amount 
of background contamination from 1.32 % to 1.1%. This directly increases the 
fit value of the lifetime by +0.2%. 

An additional systematic error is due to uncertainties in the background scale 
and Ik dependence in the three and four photon-cluster samples. Comparing 
bin-by-bin data and MC Ik distributions for background with the combined 
three and four-cluster event samples after MC signal subtraction (figure 7), 
we find an agreement at the ~ 2% level, with a small linear dependence on Ik- 
Although the agreement is at the ~ 2% level, we have conservatively taken the 
uncertainty in the overall background scale to be ±10%. The uncertainty in 
the background scale produces a systematic uncertainty in the lifetime value of 
±0.2%. Correction for the background slope changes the fit result by +0.15% 
with an uncertainty of ± 0.06% due to the statistical precision of the slope 
value. 




Fig. 7. Decay length distribution for background with three and four clusters. 
Crosses are data, solid histogram is Monte Carlo. 

Uncertainties in the DC momentum scale and the absolute EMC time scale 
enter directly in the proper time evaluation t* = Ik/ (Pjc) and give systematic 
errors respectively of ±0.1% [6] and ±7 x 10""^ (figure 1, right). 
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The fit stability has been checked by changing the lower limit of the time 
interval used in the fit between 6 ns and 12 ns and the upper between 21 
ns and 28 ns, independently. No change in is found within its statistical 
error. We have also checked the fit stability vs polar angle dividing the fiducial 
volume in two regions containing the same number of events. Specifically we 
chose events with 0.342 < | cos 6*1 < 0.766 and | cos6'| < 0.342. The values of 
tl from the two zones are consistent to within the statistical accuracy. 

In table 3 we summarize the corrections to be applied to the lifetime fit result 
and the corresponding systematic uncertainties. 



source 


correction 


systematic uncertainty 


tagging efficiency 


-0.6% 


0.25% 


acceptance 


bin by bin 


0.3% 


selection efficiency 


bin by bin 


5 • 10-5 


vertex resolution 




0.1% 


background subtraction 


bin by bin, +0.2% 


0.2% 


background shape 


+0.15% 


0.06% 


nuclear interactions 


+0.33% 


0.16% 


momentum scale 




0.1% 


time scale 




0.07% 


total 


+0.1% 


0.49% 



Table 3 

Summary of corrections and systematic uncertainties. 

The corrections add to +0.1% and the central value of the fit is moved ac- 
cordingly. The systematic error of 0.49% is at present dominated by the un- 
certainty on the dependence of tagging efficiency with Ik and by background 
subtraction. The final result is: 

tk^ = (50.92 ± 0.17stat ± 0.25syst) ns = (50.92 ± 0.30) ns 

This result differs by 1.2 cr from the other direct measurement Tr-^ = (51.54 ± 
0.44) ns [1] and by -1.8 a from the PDG 2004 fit [8], tk^ = (51.8 ± 0.4) ns. 
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